Quantum Cosmology Close to the Classical Big Bang Singularity and in the 

Semiclassical Limit 
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We investigate a cosmological model whose energy content is described by a Chaplygin gas rep- 
resented by a scalar field cf> with an associated potential producing a big bang singularity such that 
for vanishing scale factor, a — > 0, one has |<?!>| — >■ oo. The classical version of the model is discussed 
in detail, however, our main interest lies in its quantization. Upon quantization of this model in 
the Schrodinger picture, we get the Wheeler-DeWitt equation which can be solved exactly in the 
two limits a — >■ and a ~^ oo, respectively. Employing the DeWitt criterium that the wave function 
should vanish at the classical singularity in order to avoid the big bang, we show that a solution to 
the Wheeler-DeWitt equation fulfilling this condition can indeed be found. In addition to DeWitt's 
initial condition at the big bang, we postulate an asymptotic condition to be imposed on the wave 
function which guarantees that the quantum wave function is strongly peaked at the classical field 
configurations in the semiclassical limit. We also investigate a universe filled with dust (describing 
baryonic and dark matter) and show that in this case there exists an exact solution to the Wheeler- 
DeWitt equation which describes the evolution of the universe during its whole history. Also for this 
model one can construct a wave packet which avoids the big bang singularity at small scale factors 
and is strongly peaked at the classical field in the far future of the universe. Finally, we discuss the 
significance of our treatment of the Wheeler-DeWitt equation and mention some open problems. 

PACS numbers: 04.60.Ds, 98.80.Qc 



I. INTRODUCTION 

The absence of singularities may be taken to be a pre- 
requisite for any fundamental theory of nature, and in 
particular we expect a quantum theory of gravity to re- 
solve the curvature singularities found in Einstein's the- 
ory of general relativity. However, up to now there exists 
no complete quantum theory of gravity. Because of that, 
the problem of quantum cosmology has been approached 
by first carrying out a symmetry reduction (by requiring 
isotropy and homogeneity) and then quantizing the re- 
sulting minisuperspace models which have only a finite 
number of degrees of freedom 1 ■'{ . One of the main 
approaches to quantum cosmology is based on a canon- 
ical quantization of general relativity. In that approach, 
there exist basically two promising candidates for a quan- 
tum cosmological theory (see e.g. the textbooks [1] and 
[^): on the one hand, minisuperspace quantization in the 
framework of the geometrodynamical approach, and on 
the other hand loop quantum cosmology (LQC hereafter) 



the expression "singularity avoidance" means. Both ap- 
proaches describe the universe via a wave function on 
configuration space which has to be the solution of con- 
straint equations. In a canonical setting, the dynam- 
ics is implemented by the Hamiltonian constraint. The 
difference between both approaches lies in the way this 
equation is quantized. LQC is the quantization of ho- 
mogenous and isotropic minisuperspaces, where one uses 
a so-called polymer representation instead of the conven- 
tional Schrodinger representation based on the full the- 
ory of loop quantum gravity (LQG hereafter) As in 
LQG, the underlying geometry in LQC is discrete and the 
scale factor operator has discrete eigenvalues. Quantum 
dynamics is governed by a discrete difference equation in 
eigenvalues of the volume operator V [ifll • In the 
Schrodinger picture, one arrives at a differential equation, 
the Wheeler-DeWitt equation |Tl[ . In the continuum 
limit for large volume, the Wheeler-DeWitt equation is 
recovered from the loop quantum cosmological difference 
equation [12 1. 
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In the geometrodynamical as well as in the LQC ap- 
proach, one of the important problems is to investigate 
whether the classical singularities can be avoided. This 
includes that for each approach one has to define what 
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In this paper, we restrict the discussion to minisuper- 
space quantization in the Schrodinger picture. In the 
Schrodinger picture, several models have been investi- 
gated regarding their ability to cure the singularity prob- 
lem which is prevalent in the classical theory of general 
relativity. In this setting, singularity avoidance is de- 
fined as either a vanishing of the wave function at the 
point of the classical singularity (DeWitt criterium) or a 
spreading of semiclassical states indicating a break-down 



2 



of semiclassical concepts in general. In the following, we 
analyze whether the big bang singularity can be avoided 
and how reliable those solutions of the Wheeler-DeWitt 
equation are. 



We first discuss at length the Chaplygin gas [l3| which 
presents one of the simplest models possessing a big bang 
singularity. The Chaplygin gas is a perfect fluid which 
satisfies the equation of state p = —A/e, where A is a 
positive constant. The Chaplygin gas has acquired some 
popularity in cosmology as candidate for unifying dark 
sectors of the universe, that is dark energy and dark mat- 



ter, introduced by Kamenshchik et al. IJ], and has been 
widely studied since then in this context (see, e.g., [15- 
24|). Furthermore, in particle physics the Chaplygin gas 
has also raised interest due to its connection with string 
theory and its supersymmetric extension 0. The 
problem of singularity avoidance has also been investi- 
gated for an anti-Chaplygin gas for which the constant 
A above is negative 27 1. 



We also investigate a pure dust model starting with a 
big bang which describes baryonic and dark matter and 
plays the role of a prototype model for our Universe. 



This paper is organized as follows. In the next Section 
we briefly review the minisuperspace model. In Sec. IIIII 
we discuss in detail the Chaplygin gas as a simple clas- 
sical model exhibiting a big bang singularity. In Sec. IIVI 
the Wheeler-DeWitt equation for this model is studied 
and the exact solutions for the quantum wave functions 
are given for the two physically interesting limiting cases. 
In Sec. IIV Al the Wheeler-DeWitt equation is solved in 
the vicinity of the classical big bang singularity, and in 
Sec. IIV B1 in the limit where the cosmic scale factor goes 
to infinity, that is for the far future of the universe. It 
is shown in Sec. |V] that one can construct in a natural 
way solutions to the Wheeler-DeWitt equation describ- 
ing wave packets which vanish at the classical big bang 
and thus fulfill DeWitt's criterium. In Sec. lVII we present 
a detailed discussion of the wave function of the Chaply- 
gin gas in the semiclassical limit using the method of sta- 
tionary phase. We explicitly show how a wave packet can 
be constructed which is for large scale factors in leading 
approximation strongly peaked exactly at the classical 
field trajectories. Sec. IVIII is devoted to the treatment 
of the dust model for which the Wheeler-DeWitt equa- 
tion is solved exactly in the whole minisuperspace. Again 
we are able to construct a wave packet which avoids the 
classical big bang. Finally, Sec. IVIIll contains a summary 
and discussion. 



II. THE MINISUPERSPACE MODEL 

In most investigations, the problem of applying quan- 
tum gravity to cosmology is simplified by a symmetry 
reduction to homogeneous and isotropic geometries. In 
this Section, the symmetry reduction is performed on 
the classical level leaving only one gravitational degree of 
freedom given by the cosmic scale factor a{t). Restrict- 
ing to a flat Friedmann-Lemaitre universe, the metric is 
given by the Robertson- Walker metric {c — 1) 



a {t)dx 



^2 



(1) 



where we have set the lapse function equal to = 1. In 
order to take the energy content of the universe into ac- 
count, we consider the universe to be filled with a perfect 
fluid with energy density e and pressure p = p{e). At the 
quantum level, we need a more fundamental description. 
Thus, in the following, we mimic the perfect fluid by a 
homogeneous scalar field (j){t) in terms of which e and p 
are given by 



p- 



(2) 



2 " 2 

Here ¥{(/)) is the potential of the scalar field (p. In Sec. IIIII 
we shall derive an explicit expression for V{(j)) by de- 
manding that the perfect fluid is a Chaplygin gas, and in 
Sec. IVIII for a pure dust model. 

The full action, being the sum of the Einstein-Hilbert 
action and the action of the minimally coupled scalar 
field, is then given by 



S[a, (pl = / dt I 



dtC : 



dt 



,(3) 



where k = SvrG, Vo — Jj^^d^x, and = [q^^q^) := 
{a,(j)). Here we have assumed that the spatial section 
A^'^ of the universe is compact possesing the comoving 
volume VqO? . Because of the restriction of the infinitely 
many degrees of freedom of superspace to only two, one 
arrives at a two-dimensional minisuperspace defined by 
the metric 

6Vo „i 



Gab = 









Voiq'r 



(4) 



in terms of which the Lagrangian density C = C{q^,q'^) 
can be read of from From £ one obtains the canonical 
momenta tta = {tti,'^2) ■= {'^a,'^4,) 



71"! 



7r2 



dC 



dC 



_6Vo 



aa. 



(5) 
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The canonical Hamiltonian H is then given by 

niq^, tta) ig-^^TT^^B + Vo {q'f V{q') , (6) 

where denotes the inverse metric. Exphcitly, one 

obtains 

(7) 

The canonical Hamiltonian H is constrained to vanish, 
and expressing ttq and tt^ in ([7|) in terms of a, a, 4>, (p 
via ([5]) the Hamiltonian constraint becomes identical to 
the Friedmann equation 




(8) 



where H := ^ denotes the Hubble parameter. In addi- 
tion, we have the equation of motion of the scalar field 
(j), i.e. the Klein-Gordon equation 



(9) 



which is automatically fulfilled if the energy density e 
obeys the equation of energy conservation 



e = -3H{e + p) 



(10) 



III. THE CLASSICAL BIG BANG MODEL 

In the following, we assume that the perfect fluid mim- 
icked by the scalar field (f> obeys the equation of state of 
a Chaplygin gas [l^ 



P = P(e) 



A 
e 



(11) 



where A is a positive constant of dimension energy den- 
sity squared. Inserting (1111) into the continuity equa- 
tion , one obtains e in terms of the cosmic scale fac- 
tor a 



e{a) 



(12) 



where B :— (eg — A) Qq, and eo resp. uq denote the en- 
ergy density resp. the scale factor at the present epoch. 
From the Friedmann equation ^ follows that eo is iden- 
tical to the critical energy density e^rit ■— which 
is completely determined by Newton's gravitational con- 
stant G and the Hubble constant Hf). Requiring that e 
stays real and positive for all values of a > 0, we ob- 
tain < A < e^. (The case A = corresponds to a 
pure dust model, whereas the case A = respectively 



_B = corresponds to a pure cosmological constant.) Ob- 
viously, the energy density e interpolates between a pure 
dust model with e — and p = at the big bang, 
i.e. for a — > 0, and a positive cosmological constant 
A := k^/A with equation of state wa := = — 1 in 
the limit a — )■ oo. Introducing the dimensionsless param- 
eter {) < a < 1^ hy A —: aep, equation (|12p takes the 
simple form 



— = \ a 



(1 



a) [ — 



(13) 



as a function of ^ = 1 -I- z, where z is the redshift. In 
Fig. [T] we show — for o; = 0.897. As a function of z the 
equation of state w{z) -, is given by 



w{z) = — - 



a + {l-a){l + zf ■ ' 
Using ([T^ one gets from the Friedmann equation ([SJ for 



the expansion rate of the universe 



a 
do 



a 
ao 



(1 



(15) 



where do — ciqHq denotes the expansion rate at the 
present epoch. In Fig. [5] we present ^ for a = 0.897. 

Combining the two basic equations ^ and pUj) implies 
quite generally the following equation 

■d^~^ie + 3p)a (16) 
6 

for the acceleration of the universe which in the case of 
the Chaplygin gas reads 

(17) 



One concludes from ([T7)) that the universe expands faster 
and faster (d > 0) for redshifts obeying z < (j^^)^ — 
1. In order to obtain an accelerated expansion of the 
universe at the present epoch, one requires a > ^ . Fig. [3] 
shows — Ijt as a function of — for a = 0.897. 

aoHg ao 

Integration of (ITSt yields for small scale factors a{t) 
(a{t) < (f )^ = {^)iao' 
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a{t) 



'a da 



(1 + 4d6) 



from which one obtains the following expansion for the 
scale factor ait) as a function of cosmic time t near to 
the big bang (t ^ 0+) 



a{t) = ao 



1 



30 1 - a 



O 



(19) 
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where ti :— | (1 — a)"' tn and tn '■= denotes the 
Hubble time. As expected, the leading term in the ex- 
pansion (|19p agrees with the time-dependence of a pure 
dust model. 

It is worthwhile to present an alternative relation for 
the scale factor valid at all times < t < oo. Let us 
rewrite (fTSt as follows 



3 ^ 



a(t) 

ad 



da 



B 



3 Ju{a) ~1 

where we have made the substitution x 



(20) 



withC/(a) := + 



We then obtain 



t 



1 



3a 4 



arcoth C/(a) -I- arctan 



1 



Uia) 



(21) 



This relation gives for the age of the universe 
1 



tn = 



3a 4 




arctan 



(ai) 



tH ■ (22) 



At large times, t — > cx), one derives from (|21|) an expo- 
nential increase of the scale factor 



a{t) = flo 



1 — a 
8a 



exp 



(23) 



in complete agreement with the expected behavior in this 
limit due to a cosmological constant A = k\/A yielding 

y 1 

y ^ = 7^- In Fig- IH we display the evolution of the 
scale factor as a function of cosmic time for a = 0.897 
belonging to an age of the universe to = tn- 

It remains to compute the cosmic evolution of the 
scalar field (p of the Chaplygin gas and the potential V{(j)) 
associated with it. It turns out to be sufficient to know 
(j) — (p{a) as a function of the scale factor a. From ([2]) 
follows = ^py'e -|- p, and using ^ = j, the Friedmann 
equation ([8]) and the equation of state (fTTj) for the Chap- 
lygin gas combined with the explicit expression (jl2p for 
e, one obtains after one integration with 0(oo) = 

1 



0(a) 



±- 



arcoth^ 



with 



e = C(a) 



1 



A 



-a" 




(24) 



(25) 



From ([M]) and (^5]) one derives for a —t' the asymptotic 
behavior 
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FIG. 1: The energy density e of the scalar field (Chaplygin 
gas) as a function of the cosmic scale factor a for = 0.897eo- 



and we conclude that 4>{a) — )■ ±cxd near to the big bang. 
At large times, a — c», the scalar field vanishes like 



O 



V a 



(27) 



In Fig. [5] we present the two branches of 0(a) for A = 
0.897eg. 

To derive the corresponding scalar field potential V, 
we use 



(28) 



and 



e(0) = VIcosh^v/SKt/)) (29) 
which follows from ([12]) , ([24]) and ([25]), and obtain 



1/(0) 



Va 



cosh 



(\/3^0) 



1 



cosh(v3K(J 



(30) 



Note that the energy density e(0), Eq. (f^ . and the 
potential T/(0), Eq. (1501) . depend only on k and A and no 
more on the integration constant B reflecting the equa- 
tion of state (|lip as it should. The potential V{<j)) is 
shown for > in Fig. i] for A = 0.8976^. 



IV. THE QUANTUM BIG BANG MODEL 

In this Section, we shall describe the quantization of 
the classical scenario discussed above. Quantization is 
carried out in the canonical approach in the Schrodinger 
picture. Implementing the Hamiltonian constraint via 
Dirac's quantization program of constrained systems, one 
arrives at the Wheeler-DeWitt equation [ll|Tl[ depending 
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FIG. 3: Evolution of the cosmic acceleration a as a function 
of the scale factor o for ^ = 0.897eo. 



FIG. 6: The scalar field potential V corresponding to the 
Chaplygin gas depicted over ifi for the branch (p > and 



Vbb '■= exp(3\/2v') for ^ ^ oo. 



= I ^ 1.0 
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behavior at the classical singularity and in the semiclas- 
sical limit. A detailed discussion of the influence of the 
boundary conditions on the wave function is presented in 
Sec.lVl 



Solution to the Wheeler-DeWitt equation in 
the vicinity of the big bang 



FIG. 4: Evolution of the scale factor a{t) as a function of 
cosmic time t for A — 0.897en- 



on the configuration space variables (a, 0) (respectively 
(x, if) to be defined below). From the solution obtained in 
this Section, we shall retrieve information concerning the 
above Chaplygin gas model with regard to its quantum 



Now we analyze the quantum version of the classical 
model discussed in Sec. IIIII The wave function satisfies 
the Wheeler-DeWitt equation 



2 



K 



^ (x, <p) = 



where we have introduced the dimensionless variables 
X := In^ and ip := yf\<\> and used the Laplace-Beltrami 
factor ordering. The potential Vi^p) is given by Vi^p) '■= 



6 



V (^1'/') > see equation pi)) . 

In this Section, we are interested in the vicinity of the 
big bang, where one has — >■ c» and thus we can ap- 
proximate the potential by an exponential function, see 
Fig. ini In the following, we study the case > 0. Anal- 
ogous results hold for the branch <p < 0. Thus we obtain 
near to the big bang the Wheeler-DeWitt equation : 



Here the functions 



92 



a- 



2 ; *bb (a;, (p) 



+ ^^VlaleiC£,\2,(2x + V2ip\\ ^-bb (a;, V') = 0. (32) 
ZK L V /J 

Here the index bb denotes the approximation to the wave 
function near to the big bang (bb). We introduce new 
coordinates u = u{x,ip),v — v{x,(p) (see also [27]) in 
order to reduce Eq. ([5^ to 



fi2 

where 



g2 Q2 
^-^ - TTT ) 'J'bb {u, v) + *bb {u, v) 







'.{x,Lp) := 



A r 



pe 



and 



A := 



(33) 

(34) 
(35) 



2Ai 



'3k 



= Vna: 



(36) 



The solution to psp can be found from the ansatz 

*bb {u, v) = /(71W -I- 72w) , (37) 

where 71 and 72 do not depend on {u,v). Inserting (j37p 
into the Wheeler-DeWitt equation ([55]) leads to 



(38) 



where the prime denotes a derivative with respect to 
the whole argument of /. Setting ft^ (7J — 7I) = 1, we 
get with 72 := ±|, > 0, 71 = ^ and uj — uj{k) 
+y/k'^ + 1 the simple differential equation 



/" + / = 



(39) 



which has the two solutions (C^^ are arbitrary complex 
coefficients) 



*bb(", v) = C± sin ( % ) + C± cos ( ^ ) . (40) 



h 



S± ~ S±{k; u, v) := uj{k)u ± kv 



(41) 



denote the classical actions satisfying the Hamilton- 
Jacobi equation corresponding to the Wheeler-deWitt 
equation ([33]). Although the solutions (l40l) solve the 
Wheeler-DeWitt equation for every complex value of k 
(with u){k) being then complex, too), we have chosen the 
parameter k to be real and positive, since the physical 
boundary conditions to be discussed in the next Sections 
necessarily lead to the interpretation of fc as a real wave 
number. Furthermore, we note that the two other ad- 
missible classical actions obtained by S± —S± do not 
lead to new linearly independent solutions of (I33|) . Thus 
we conclude that the complete solution of the Wheeler- 
DeWitt equation is given by ^f+j^ + -^^^y^. In Sec |V] 
we shall show that the coefficients Cj^j determined 
by the physically correct boundary conditions . 

Here we want to point out that the exponential approx- 
imation to the exact Chaplygin gas potential (|30t does 
not interpolate between matter and dark energy, contrary 
to the Chaplygin gas. In fact, it mimics a pure dust com- 
ponent with p = respectively w — (see Sec. lVIip . This 
means, due to the fact that the pressure of the Chaplygin 
gas behaves like p ~ — -^a^ + . . . near to the big bang, 

that the exponential potential Vbb describes the Chaply- 
gin gas precisely at the big bang, a = 0. Including higher 
order terms in the expansion of the Chaplygin gas poten- 
tial in the vicinity of the big bang would be interesting, 
but would lead to a more complicated Wheeler-DeWitt 
equation, whose exact solution has not yet been found. 



B. Solution to the Wheeler-DeWitt equation for 
large values of the scale factor 

In the case of large values of the scale factor, a 00, 
respectively x ^ 00, the classical field ip approaches zero, 
see Eq. (P7)) , and therefore we can approximate the scalar 
potential V{(p) in the Wheeler-DeWitt equation (|3T|) by 
its leading term ^/A, since V{ip) = V^(l + 0{ip^)) for 
93 0. Thus we study in the following the Wheeler- 
DeWitt equation 



-H9A2e6-*„o {x.p) 



, 



(42) 



where the index 00 denotes the case a — 00. To solve 
this equation, we make the separation ansatz 



(43) 
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and arrive at 



Coo (2^) 



respectively 



K 



$00 M 



(44) 



(45) 



Thus the solution for the matter-dependent part of the 
wave function "^^a is given by 




In order to solve the gravitational part of Eq. (H^ 



(46) 



FIG. 7: The classical trajectory <f^\ (positive branch) as a 
function of the coordinates u and v. The dashed lines v = ±u 
indicate the boundaries of the physically allowed sector. 



.6a; 



K 



Coo(a:)=0, (47) 



AVOIDING THE BIG BANG SINGULARITY 



we make the ansatz 



Cooix) = F{z), 



(48) 



where z := f3e^^ . Inserting Eq. (I15|) into (|47p . one gets 
the Bessel differential equation 28] 

z^F"{z) + zF'{z) + (z^ _ y-^) F{z) = 0, (49) 

if (3 is chosen to be equal to (see Eq. ([M]) ) and is 
given by i^. Thus we get the following solution for the 
gravitational part: 



^3x 



, (50) 



where J„{z) is the Bessel function and Yi,{z) the Neu- 
mann function. We can rewrite the above wave function 
in terms of the Hankel functions hI^"^ (z) 



Coo{x)^ViH'-£ I -e' 

3/1 \ ll 



A 



3a; 



2?2i?fi ( Je3-). (51) 



A 



The general solution to the Wheeler-DeWitt equa- 
tion (|42)) is then given by: 



V^H^l \ -e^ 



(52) 



with the arbitrary constants Ai , T^i and K which will be 
determined in Sec. I VII 



In Sec. IIV Al we have derived the general solution (^0)) 
to the Wheeler-DeWitt equation ([32)) respectively ([33)) 
being an approximation valid near to the classical big 
bang to the exact Wheeler-DeWitt equation ([3T)) . This 
solution still depends on the four coefficients 
the constant k. It is the purpose of this Section to show 
that a physically sensible and mathematically consistent 
choice of the boundary conditions close to the big bang 
and in the far future exists, that avoids the classical big 
bang singularity and is strongly peaked exactly at the 
classical field solution in the semiclassical limit. 

The solution ([10)) has been found in terms of the aux- 
iliary variables u and v defined in Eqs. p4p and ()35l) . 
In order to obtain their physically allowed range, one 
must go back to minisuperspace with the coordinates 
(a, iy9) € ]R> X Il>. (Remember, that we restrict the 
discussion to the positive branch of the scalar field ip). 
It then follows from ([M)) and (1551) that the physically 
allowed domain of u and v is given by the sector S := 
{{u,v) G Il> X M, \v\ < u} whose boundary is indicated 
by the dashed lines in Fig. [7| We also show in Fig. [7] the 
classical trajectory ipd := (pciiaiu,v)) = y/^(f){a{u,v)) 
given in Eq. ()24p (with the upper sign) and shown in 
Fig. [5] as a function of the scale factor a. 

As is well-known, there exists no genuine time variable 
in the geometrodynamical approach. One can use, how- 
ever, either the cosmic scale factor a or the scalar field 
if as an "effective time" variable (emergent time). In 
the following, we consider the scale factor a (respectively 
X := ln(^)) to play the role of a "clock" for the quan- 
tum evolution. The deep quantum regime corresponding 
to the classical big bang is then defined as the hypersur- 
face a = 0, and thus the quantum evolution requires the 
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definition of an initial condition imposed on the quan- 
tum wave function on this hypersurface (a = 0). Here, 
we choose De Witt 's initial condition on the hypersurface 
a = in the minisuperspace (a, ip) G R> x Hy 



Van 5* (a, ip) 

a— J-O 



0, 



(53) 



which should hold for any given value of the scalar field 
(fi E R>. Note that we do not require that ip has to be 
evaluated on the classical trajectory cpciia)- However, 
the condition (|53p does not completely determine the 
wave function, but rather has to be supplemented by an 
asymptotic boundary condition controlling the behavior 
of 4* in the far future, a — >■ oo, as discussed below. 

We are now in the position to give the initial condition 
on the wave function ^>{u,v), Eq. (|40p . considered as a 
function of u and v, by translating the initial condition 
([55)1 into the {u, w)-plane. From ([Ml) and ([55)) one obtains 

u = O (^{-^y) and V ^ -u + O ((^)'") ^ for 
a ^ 0, (f fixed, and thus the initial condition becomes 



lim ^'(w, —u) = 0. 



(54) 



Using the solutions ^'bb(^i Eq. ([40]) . the initial condi- 
tion ([54)) gives C2 =0, and we obtain 

*bb(u,«) = C+sin(^^^ +Crsin(^^^ (55) 

as the solution of the Wheeler-DeWitt equation ([55)) sat- 
isfying DeWitt's criterium ([55)) for avoiding the classical 
big bang singularity. 

Let us remark that the wave function ([55)) satisfies the 
Wheeler-DeWitt equation near to the classical big bang 
together with the initial condition (|53p for any complex 
value of k. This is important since the value of k respec- 
tively its admissible values (discrete or possibly continu- 
ous) cannot be determined from a solution whose validity 
is restricted to small values of the scale factor a. To com- 
pletely fix the spectrum of k (and thus of w) requires the 
knowledge of the wave function in the whole minisuper- 
space which in turn can only be determined (at least in 
principle) if the allowed space of functions equipped with 
an appropriate scalar product is given. While we shall 
present in Section I VII the quantum wave function satis- 
fying the approximate Wheeler-DeWitt equation valid in 
the classical regime, i.e. for a ^ 00, and thus for late 
times, the complete solution for the Chaplygin gas valid 
in the whole minisuperspace is not known. 

As already mentioned, the initial condition (I53p has to 
be supplemented by a second boundary condition. Here, 
we propose the following semiclassical boundary condi- 
tion which selects from a wave packet based on the wave 
function ()55|) one which in the limit of ?i — > is strongly 



peaked at the classical field configuration (^s^ correspond- 
ing to the positive branch of the scalar field. (Remember 
that in the Wheeler-DeWitt equation ([5^ we have used 
the approximation Vbb to the scalar potential valid for 
ip > 0). Since the Wheeler-DeWitt equation is linear, we 
can form from ()55|) the most general wave packet using 
the superposition principle and obtain 



f°° f S 

*bb(M,w) = / Ci+(fc)sin( - 



.00 

/ (fc) sin 
Jo 




(56) 



(Here we have assumed that the spectrum of k is contin- 
uous and given by fc > 0). Note that Eq. ([55)) displays 
explicitly the dependence on ?i, since the classical actions 
S± are independent of h. The crucial point is that the 
wave packet ([56)) still satisfies the Wheeler-DeWitt equa- 
tion ([55)) with the initial condition ([51)) for any given 
complex "wave number amplitude" C^^^i^)- 
Rewriting ([56)) as follows 



*bb(w, w) 



C+(fc) 

cr(fc) 



e " — e 1 



■ s_ 



dk 



dk 



(57) 



we observe that ([57)) has the standard form of a WKB- 
wave function which often is used as an ansatz to obtain 
an approximate solution for quantum wave functions in 
the semiclassical limit when h goes to zero. Thus the 
wave function (I57|) has a simple physical interpretation 
in terms of the classical actions S± which justifies our 
assumption of k being real. It is important, however, to 
observe that ([57)) has been derived from our exact solu- 
tion ([55)) to the Wheeler-DeWitt equation ([55)) satisfying 
the crucial initial condition ()54p without assuming the 
WKB-approximation. We would like to mention that a 
related wave function has been introduced in 27] in the 
discussion of the anti- Chaplygin gas, in that case, how- 
ever, as an WKB- ansatz which does not satisfy the initial 
condition ([5^. 

In order to investigate the limit h —>■ oi the wave 
packet ([57)) , one uses the method of stationary phase (see 
the detailed discussions in Sections IVl) and IVlI)) and finds 
that (under suitable conditions on the functions Cf{k)) 
the first integral is in the semiclassical limit strongly 
peaked at the positive branch of the classical field and 
the second one at the negative branch. Employing now 
the asymptotic boundary condition that the wave packet 
should reproduce only the classical result corresponding 
to the positive branch, we conclude Ci{k) = 0, and thus 
the exact wave packet solution to the Wheeler-DeWitt 
equation ([55]) satisfying both boundary conditions takes 
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the simple form 



*bb(u,w) 



dk 



(58) 



where we have redefined the ampliude C{k). 

The details of the calculation using the method of sta- 
tionary phase will not be discussed at this point, since 
the derivation is identical to the one to be presented in 
Sec. lVIIl where we shall study a pure dust model instead 
of the Chaplygin gas. Since the Wheeler-DeWitt equa- 
tion can be exactly solved for all values of x in the case 
that the matter content of the universe consists only of 
a dust component, we are then able to discuss the exact 
wave function in the whole kinematically allowed range. 

In conclusion, we have demonstrated that the Wheeler- 
DeWitt equation possesses near to the classical singular- 
ity a unique solution corresponding to the quantum wave 
function (j58|) which avoids the big bang singularity in the 
sense of DeWitt's criterium ([55)1 respectively ([5l)) and is 
strongly peaked at the classical field Lp'^^ {x) in the semi- 
classical limit. 



VI. QUANTUM COSMOLOGY OF THE 
CHAPLYGIN GAS IN THE SEMICLASSICAL 
LIMIT 

In equation ([5^ we have already given the general so- 
lution ^'oo(a:, (fi) to the Wheeler-DeWitt equation for the 
Chaplygin gas in the small-field approximation. To ob- 
tain from this expression the leading term in the semi- 
classical limit, we have to investigate the behavior of the 
Hankel functions H':]f ^ (-^ e*^^) for ft —> and x, K fixed. 



This limit corrsponds to the Debye type asymptotics of 
{z), where both /i := J| and z := e'^^ are large 
and real positive numbers such that ^ — Jje"'^^ is fixed 
as /i, z — )■ oo. Using the well-known asymptotic formula 
in this limit (see [28|), we obtain for ft 



^oo{x,ip) 



I 



(59) 



Bie 



(l + 0(ft)). 



Here the functions 

St = St{K; X, ip) ■.=^^K^ + m^x) 

K J K \ (60) 
— — arsmh — — ± Kip 
3 \m{x) J 

denote the two classical actions satisfying the Hamilton- 
Jacobi equation corresponding to the Wheeler-DeWitt 



equation p2)) . Furthermore, we have defined 

m{x) := SAe^'^ (61) 

and have introduced the new (iiT-dependent) coef- 
ficients 61,2 := ^e±^=F'i^i,2X'i,2 and := 
v|Ie±^T'i^2,i2?i,2- Since the Wheeler-deWitt equa- 
tion is linear, the leading term of the most general 
wave packet is given in the semiclassical limit by 

^ , , r dK 



00 ^/2^{K^ + m'^{x)Y 
Bi{K)ci'^^ + B2{K)e-Ti^^ 

+B3iK)e^^^ + B4iK)e-i^^] (1 + 0{h)) , 

(62) 

where Bi{K) now denote general complex "wave ampli- 
tudes" . We shall now show, using the method of station- 
ary phase, that the wave packet (p^ is for ft —j' strongly 
peaked at the classical field trajectory ipciix) if the am- 
plitudes Bi{K) are strongly peaked at the "classical wave 
number" Kd to be determined below. 

In a first step let us assume that the amplitudes 
Bi{K) are smooth functions with compact support and 
that the actions 5^ possess a finite number of isolated 
non-degenerate stationary points = K^ixjip) S 

supp Bi{K), 



dK 

9^2 



{K^{x,(p);x,(p) = , 
{K^ix,ipy,x,p) ^0 . 



(63) 



We then obtain for the first term in the integral (|62l) . 
denoted by '^^ooixjip), in the semiclassical limit, ft 

M 



{K+{x,(p);x,(p) 



{K+{x,p))^ + m^x) 



X exp<^ ^5+ {K+{x,p);x,ip) 



(64) 



-sgn 



dK^ 



{K+{x,ip);x,ip) 



X (l + 0(ft)) 



For the remaining three terms of (|62p we obtain similar 
expressions with the obvious substitutions made in (j64l) . 
With the classical actions given in we obtain 



dSt 1 . ^ fK\ ^ 



(65) 
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and 



and thus the only stationary points are 



K^{x, if) — ±m{x) sinh(3(/3) 



leading to 



1 



3rn(a;) cosh(3iy9) 



(66) 



(67) 



(68) 



In the last step let us assume that the functions 
Bi{K) are strongly peaked at the classical wave num- 
ber K = Kci which in turn implies that the amplitudes 
Bi [K^{x,ip)) are strongly peaked at K^{x,ip) = Kd- 
From this condition and (|57)) we conclude that e.g. 
the wave packet 'i!^{x,ip) derived in equation (1641) is 
strongly peaked at those (a;, (y9)-values which are given by 



vtA^) — ^arsinh 



m{x) 



= — arsinh 
3 



3A 



(69) 

Since ^'^■'(a;, 1^9) is a solution to the Wheeler-DeWitt 
equation for the Chaplygin gas in the small field approxi- 
mation respectively in the limit x — lii(^) 3> 1, we infer 
from ([5^ for the semiclassical field f'^^ix) the asymp- 
totic behavior 



O 



V a 



a > ao, (70) 



in complete agreement with the exact behavior (|27p in 
this limit (for the positive branch of the scalar field) if 
the classical wave number Kd takes the value 



6B 



Vo 



(71) 



Considering in a similar way the contributions to the 
wave packet "^ocix^ip) involving the action S*^, we con- 
clude - using the same value ([71]) for Kd - that they are 
strongly peaked at K~(x,ip) — Kd and thus at those 
(x, (/3)-values which are given by 



'^sci(a;) 



-arsinh | ^^^e 
3 V 3A 



(72) 



again in complete agreement with the exact asymptotic 
behavior (|27l) for the negative branch of the scalar field. 

It is interesting to compare the semiclassical results 
(ISgi) and dH]) with the exact formula ([M]), valid at all 



a > 0, which can be rewritten in terms of Kd and reads 
(for the positive (-I-) respectively the negative (-) branch) 



1 



3^2 



arsinh 



(73) 



While the leading terms in (^^j and (yS^ agree, it is seen 
from ([73]) that the higher order terms in (^) are differ- 
ent. 

We have thus demonstrated that there exists a quan- 
tum wave packet derived from the Wheeler-DeWitt equa- 
tion which in the semiclassical limit is strongly peaked 
exactly at the classical field configurations. 



VII. QUANTUM COSMOLOGY OF A PURE 
DUST MODEL 

In the previous Sections, we have investigated quan- 
tum cosmology for a universe whose matter content is 
assumed to consist of an ideal fluid satisfying the equa- 
tion of state of a Chaplygin gas. The latter is described 
by a homogeneous scalar field whereby the dynamics is 
completely fixed by the scalar potential (15(11) . The Chap- 
lygin gas interpolates between pure dust at the big bang 
and dark energy described by a positive cosmological con- 
stant at late cosmic times. 

In this Section, we study a pure dust (d) model sat- 
isfying for all scale factors a > the equation of state 
Pd = 0, i.e. a pressureless fiuid. We know that the en- 
ergy content of our Universe is dominated during a long 
period of its evolution by dust consisting of baryonic and 
dark matter. The equation of energy conservation (jlOp 
determines completely the a- respectively x-dependence 
of the energy density ed for dust which reads 

e, = eo{^y=eoe-'-. (74) 

Here eo is given, as in the previous Sections, in terms 
of the Hubble constant Hq since the Friedmann equa- 
tions ^ leads to eo = Ccrit := Furthermore, it 
follows Vd = ^ from the two equations in 1^ which then 
yields for the scalar potential 



-3x 



(75) 



where we have introduced again the dimensionless field 

2 



In addition, we obtain from ([2]) ^^^^ 

d = y/^edciox, and thus 



r ed and from 



y dx J 

scribing pure dust comes in two branches (f^ which are 
for all X € R exactly given by 



i. It then follows that the scalar field de- 



(76) 
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Here ipf (0) denote arbitrary initial values at the present 
epoch which we set in the following equal to zero without 
loss of generality. Combinig ([75]) and ([75|) . we conclude 
that the dust model is governed by an exponential po- 
tential which reads for the two branches ip^ 



(77) 



We observe that ([77|) agrees (for the branch (p^)with the 
potential Vbb of Section IIV Al used there as an approxi- 
mation to the potential (pO)) for the Chaplygin gas close 
to the big bang and shown in Fig. [5] as the dashed curve. 
(To obtain this agreement, we have to replace the param- 
eter A of the Chaplygin gas by 4eQ.) There are, however, 
important differences between the potentials Vbb and : 
i) whereas Vbb is only an approximation to the exact po- 
tential (|30|) in the large-field limit {(p — > oo), the potential 
V^"*" is exact for dust, ii) in the case of the Chaplygin gas, 
the infinite future of the universe {x — >■ oo) is reached 
for (/J — ^ 0, while for the dust component one requires 
93 J — > ^00 for a; — > 00. This is obvious by rewriting the 
energy density (j74p in terms of ip^ 



eoe 



(78) 



Thus the physically relevant domain in minisuperspace is 
(x, ip^) G R X K, for dust in contrast to (x, p^) G R x R;^ 
for the Chaplygin gas. 

From ([3T|) and ([77|) we conclude that the wave function 
^'d {x, ip) for a universe filled with dust obeys the exact 
Wheeler-DeWitt equation 



^2 (d^ _ 

\dx^ dpj\) 



{x.ipd) 



^^2g3(2.+y2^,)^^ (X,^d)=0 



where the constant Ad is defined by 



2V2H0 



(79) 



(80) 



(Equation (j79|) is written for the branch p'^ using the 
potential V^^ of Eq. ([77]) . and the superscript -I- is sup- 
pressed to simplify the notation.) A comparison with the 
Wheeler-DeWitt equation (l32l) reveals that (|79l) can be 
recast into the form 

^'(^|^-^)*d(w,«) + *d(w,«) = (81) 

using the new coordinates (u, v) which are defined as in 
equations (p4| - ([36)) . but replacing A by Ad given in ([SO]). 
Although the range of (a;, (^) is R x R in ([75)1 instead of 
(R X R> ) in ([5^ , it is obvious that the physically allowed 



domain of (u, v) for dust is again given by the sector 
S := {{u, u) e R> X R, \v\ < u} indicated in Fig. [T] 

It follows that the most general exact solution to the 
Wheeler-DeWitt equation (fSTj) for dust satisfying at the 
classical big bang singularity DeWitt's criterium (|53|) re- 
spectively ([54|) is given by (see equation ([57])) 



Vfd («,!')= / C+(fc) 



Cd(fc) 



dk 



e' R 



(82) 



Here the classical actions S± are defined in (|41l) with 
fc > and uj{k) +^/k'^ + 1. C^{k) denote arbitrary 
complex " wave number amplitudes" . Although the wave 
packet (|82|) has the form of a WKB- ansatz, it is important 
to remark that ((5^ is in the case of dust the exact solu- 
tion to the Wheeler-DeWitt equation ([8T|) . Furthermore, 
in contrast to the solution (157)) for the Chaplygin gas, 
which is only an approximation close to the big bang, 
the wave packet ([5^ is for a universe filled with dust 
valid for all cosmic scale factors a > 0. Since 5'd(w, w) 
in (15^ vanishes for a — > 0, that is if one approaches the 
classical big bang singularity, the latter is avoided in the 
sense of DeWitt's criterium. It remains to investigate the 
wave packet ((82|) in the semiclassical limit in order to see, 
whether the amplitudes Cf{k) can be chosen in such a 
way that the wave packet is for h ^ strongly peaked 
at the positive branch p'^{x) of the scalar field. 

Evaluating the wave packet (|82l) in the semiclassical 
limit using the method of stationary phase (assuming 
C^{k) to be smooth functions with compact support), 
one obtains up to corrections of order h that the two 
integrals in ([82]) are proportional to Cj"(fc+(u, w)) respec- 
tively C^(fc:^(M, u)) (see the discussion in Section |VI|. 
where k^ = k^{u,v) e suppC^(A;) denote the stationary 
points determined by 



dS"^ {k^{u, v); u, v) 
dk 

d^S^ {k^iu,vy,u,v) 



k^{u,v) 



(fct(",«)) +1 
u 

{k^{u,v)f + 1 



(83) 



^0 



The only stationary points are then given by 



k^ {u,v) 



(84) 



^ 0, keeping in mind that 



where k^ corresponds to v ^ 

k > must hold. Furthermore, the second derivative in 

_ 3 

^ 2 



83t takes the value u 



1 



which is nonzero in 



the interior of the physical region, i.e. for < < u. 
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Let us now assume that both amphtudes C^{k) are 
strongly peaked at the classical wave number 

feci 1 . (85) 

It then follows that the wave packet \I'd(it, u) is in the 
semiclassical limit strongly peaked at :— which 
in turn implies that it is strongly peaked precisely at the 
classical field configurations ip^ix) = T"^, in agreement 
with the exact solutions (|76|) (for (0) — 0). The last 
statement is easily shown to be true, e.g. for (p'^{x), by 
calculating m+ := u{x,ip'^) = {a + p)^e~^ and w+ := 
vix^ip-^) = -{(7 - p)^e^'' and noting that ^ 

which implies the correct relation = — 

v2 

It is now obvious that one can fulfill the remain- 
ing boundary condition which guarantees that the wave 
packet ([82|) is in the semiclassical limit only peaked at 
the positive branch of the scalar field (corresponding to 
the choice of the potential V^{ip^) considered in the 
Wheeler-DeWitt equation ((79)) ) by imposing the condi- 
tion C^(fc) = 0. 

It has thus been demonstrated that in the case of a uni- 
verse filled with dust there exists a quantum wave packet 
derived from the exact Wheeler-DeWitt equation which 
avoids the classical big bang singularity and is strongly 
peaked for fi. — ^ precisely at the exact classical field ip'^ , 
valid at all cosmic scale factors a > 0. 

VIII. SUMMARY AND DISCUSSION 

In this paper, we have studied a Friedmann-Lemaitre 
universe whose energy /matter content is described by a 
scalar field (j) obeying either the equation of state of a 
Chaplygin gas or that of pure dust. In Sec. |TT]the gen- 
eral symmetry reduction to homogeneous and isotropic 
geometries has been performed in the frame work of clas- 
sical cosmology starting from the Einstein-Hilbert action 
and the action of a scalar field minimally coupled to grav- 
ity and governed by an arbitrary potential V{(j)). The 
canonical Hamiltonian of the two-dimensional minisuper- 
space (a, 0) obtained this way is constrained to vanish 
and leads to the standard Friedmann equation for the 
cosmic scale factor a{t) and to the Klein-Gordon equa- 
tion for <p{t). Specializing to the case of a Chaplygin 
gas completely determines the scalar potential V{(j)), see 
equation (j30l) . The classical time evolution of a{t) and 
(j){t) has been investigated in detail in Sec. IIIIl and it is 
seen that the universe starts in this model with a classical 
big bang as in a pure dust model and undergoes at late 
times an accelerated expansion dominated by a dark en- 
ery component as in models with a positive cosmological 
constant. The results of this model have been illustrated 



in Figs. [T][7| for a particular value of the parameter A (see 
equation pip ), being the only free parameter character- 
izing the Chaplygin gas, which is chosen in such a way 
that the age of the universe is equal to the Hubble time. 

In Sec. IIVI the quantization of this model has been 
carried out in the framework of geometrodynamics lead- 
ing to the Wheeler-DeWitt equation for the "wave func- 
tion of the universe" 5'. For a general scalar field, the 
Wheeler-DeWitt equation in the two-dimensional minisu- 
perspace choosing the Laplace-Beltrami factor ordering 
has been given in equation (pij) using the dimensionless 
variables x := In and ip := y^4>- Since the ex- 

act solution to the Wheeler-DeWitt equation valid in the 
whole minisuperspace is not known to us, we have studied 
in Sec. IIV Al the approximation to the Wheeler-DeWitt 
equation which holds close to the classical big bang singu- 
larity (corresponding to \(p\ — ?> oo), and in Sec. IIV B] the 
alternative approximation valid for large scale factors re- 
spectively small scalar fields ip 0. In both limits, the 
complete solutions to the Wheeler-DeWitt equation have 
been found, see (gH]) respectively equation ([5^. These 
solutions depend on arbitrary coefficients which have to 
be fixed by boundary conditions as one is familiar from 
standard quantum mechanics. 

In Sec. E] we have shown that the classical big bang 
singularity can be avoided if one imposes DeWitt's ini- 
tial condition (1531) ascertaining that the wave function 
vanishes on the hypersurface a = for any given scalar 
field ip. Under this condition, the remaining freedom in 
the wave function is already drastically reduced as shown 
in equation (j55p . To reduce this freedom even further, we 
have proposed to supplement DeWitt's boundary condi- 
tion at the big bang by a second boundary condition to 
be imposed in the semiclassical limit in which Planck's 
constant approaches zero. Since in this limit the stan- 
dard classical cosmology should be recovered, we demand 
that the wave packet constructed from the solution ((55)) 
should be strongly peaked (with a suitable choice of the 
wave number amplitude) at the correct classical solution 
of the scalar field ip. In this way we then obtained the fi- 
nal solution ((58)) of the Wheeler-DeWitt equation which 
has the form of a WKB-wave function even though it 
solves the Wheeler-DeWitt equation exactly. The fact 
that the wave packet ((55)) depends only on the classical 
action S+ (corresponding to the classical outgoing field) 
automatically guarantees that our second condition is ful- 
filled in the semiclassical limit as we have shown by the 
method of stationary phase. 

In Sec. I VII we have explained in detail how the method 
of stationary phase is applied in the semiclassical limit 
to the wave function ((52)) which represents the complete 
solution to the Wheeler-DeWitt equation ((^^ in the case 
of large values of the scale factor, that is for the universe 
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in the far future. Employing the Debye type asymptotics 
of the Hankel functions, we have constructed the gen- 
eral wave packet ([5^ which gives the dominant term in 
an asymptotic expansion for /i — > 0. Under the assump- 
tion that the wave number amplitudes Bi{K) are smooth 
functions with compact support and are strongly peaked 
at the classical wave number Kc\ (given in Eq. (|7ip ). 
we have shown that the wave packet is strongly peaked 
at the semiclassical field trajectories ^f^ii^) respec- 
tively (17^ . The main point to observe is that these semi- 
classical fields possess for a ^ oq in leading approxima- 
tion exactly the same behavior as the exact classical fields 
'Pel i^) given in ([75| . 



Finally, in Sec. IVIII we have studied the quantum 
cosmology of a pure dust model which is governed by 
the exponential potentials (|77p corresponding to the two 
branches of the classical field. In this case the Wheeler- 
DeWitt equation ([7^ holds in the whole minisuperspace 
without any restriction. Imposing DeWitt's initial con- 
dition ([55]) , we have derived the general wave packet ([5^ 
satisfying the Wheeler-DeWitt equation and have thus 
found the exact quantum wave packet for a universe filled 
with dust which avoids the classical big bang singularity. 
Taking into account the second boundary condition in the 
semiclassical limit, we have shown that from ((82|) one can 
construct a wave packet which is strongly peaked e.g. at 
the exact field ip'^ (x) corresponding to the dust potential 
V+ in dZll). 



We would like to remark that the various solutions 
to the Wheeler-DeWitt equation have been obtained 
under the assumption that this equation can be treated 
as in ordinary quantum mechanics. The question 
whether the Wheeler-DeWitt equation corresponds to a 
self adjoint operator in a Hilbert space equipped with 
a scalar product leading to a consistent probability 
interpretation has not been dicussed in this paper. To 
our knowledge, there exists up to now no complete math- 
ematical treatment concerning these problems in the 
framework of geometrodynamics (see however 3^, 31|). 
It is therefore of great importance that there exists an 
alternative approach toward the quantization of gravity, 
namely LQG respectively LQC if applied to cosmology. 
In this theory the avoidance of singularities could be 
accomphshed according to 0-0, S|- Moreover, the LQC 
ansatz can be used to make physical predictions in the 
vicinity of the big bang. As shown in j32| the situation 
is improved with respect to the singularity avoidance 
in scalar field models such as the (anti-)Chaplygin gas, 
inflationary, radiation dominated or quintessence models. 
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